Infrared and optical properties of pure and cobalt-doped LuNi2B2C. 
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We present optical conductivity data for Lu(Nii_a;Coa;)2B2C over a wide range of frequencies 
and temperatures for x=0 and x=0.09. Both materials show evidence of being good Drude metals 
with the infrared data in reasonable agreement with dc resistivity measurements at low frequencies. 
An absorption threshold is seen at approximately 700 cm~^. In the cobalt-doped material we see 
a superconducting gap in the conductivity spectrum with an absorption onset at 24 ± 2 cm~^ = 
3.9±0.4fcsrc suggestive of weak to moderately strong coupling. The pure material is in the clean 
limit and no gap can be seen. We discuss the data in terms of the electron-phonon interaction and 
find that it can be fit below 600 cm^^ with a plasma frequency of 3.3 eV and an electron-phonon 
coupling constant At,- = 0.33 using an o? F{ijj) spectrum fit to the resistivity. 
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I. INTRODUCTION 

The recent discovery of superconductivity at 39 K in 
MgB2 has renewed interest in intermetaUic superconduc- 
tors containing boron. There have been suggestions that 
their high transition temperatures may result from an ex- 
otic mechanism and not from a conventional s-wave BCS 
phonon process. The borocarbide family of superconduc- 
tors LNi2B2C, where L=(Y,Lu,Tm,Er,Ho,and Dy) shows 
a number of physical properties that suggest that they 
may be a testing ground for these ideas.EJ For example, 
recent measurements of heat capacity and microwave sur- 
face impedance in the vortex state suggest that law en- 
ergy quasiparticles may have d-wave dispersion.ai3 The 
critical field temperature variation has an upward curva^ 
turc that is difficult to explain within the BCS theoryBcl 
Other suggestions of abnormal behavior include reskU 
ual absorption in the gap seen by Raman spectroscopyBcl 
and the sensitivity ot|the transition temperature to non- 
magnetic impurities .Lfu On the other hand, there is a 
lot of the evidence that supports a conventional s-wave 
mechanism of superconductivity for these materials. This 
includes good agreement between transport atud—criti- 
cal field properties based oup^liashberg theoryaEj and 
s-wave-like tunneling spectraO 

Within the s-wave, electron-phonon picture, there are 
two routes to high transition temperatures: a large cou- 
pling constant A or the coupling to very high frequency 
modes with modest coupling. It was suggested early on 
by Pickett et al. that unusual electronic structure seen 
in the borocarbides at the Fermi surface would lead to 
strong coupling to certain phononsilj On the other hand 
more recent data based on specific heat by Michor et al. 
axe consistent with a A of the order of unityJlj Gonnelli 
et al. were able to fit the temperature dependence of 
the dc resistivity of YNi2B2C with a phonon spp:trum 
from inelastic neutpQn scattering of Gompf et al.lB with 
a \tr of only 0.57J13 where \tr is just A modified by a 



factor that depends of the scattering angle. □ Gompf et 
al. found evidence of soft mode behavior in several modes 
but not th&jhigh frequency mode originally singled out by 
MattheissO to couple strongly to the electronic system. 

The infrared response of superconductors can be used 
to sort out some of these issues. For example, there are 
striking differences between exotic and conventional su- 
perconductors in the effect on the infrared response of 
doping with impurities that limit the scattering time of 
the quasiparticles. Doped d-wave superconductors do not 
develop a sharp gap signature at a; = 2A as they ap- 
proach the dirty limit where the scattering rate l/r > 2 A 
as described by Mattis and Bardeen.H Instead of a region 
of perfect reflectance followed by a sharp onset of absorp- 
tion at 2 A, as seen in conventional materials, d-wave ma- 
treials show a low-frequency Drude-like absorption fram. 
defects surrounded by regions of normal material. E£lE2l 
For this test it is necessary to measure the reflectance 
accurately in the gap region of doped samples in the re- 
gion of huj = 2A. 

It has been suggestecD that the sensitivity of Tc to non- 
magnetic defects in the borocarbides is due to changes in 
N{0), the density of states at the Fermi level, that are 
related to a peak in the density of states near theJ'ermi 
energy Ep seen in band structure calculations £30 One 
effect is the valence effect, where, in the rigid band pic- 
ture, doping with holes or electrons will sweep Ep across 
this peak and change N{0). Another is where doping 
with non-magnetic impurities broadens the peak which 
leads to the reduction of N{0) which in turn reduces Tc. 
One of the advantages of infrared spectroscopy over dc 
transport measurements is that it allows a separate de- 
termination of the scattering rate l/r and the free car- 
rier plasma frequency cOp which is related to N(0). Any 
dramatic changes to the plasma frequency with doping 
would be seen as a change in the Drude spectral weight 
with doping. There will also be changes to the slopes of 
the resistivity curves with temperature which, in accord 
with Matthiessen's rule, should not change with doping 
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if it affects tlie carrier life time but not the density of 
states. 

Infrared spectroscopy can also be used to map out the 
spectrum of excitations that couple to the electrons. It 
was shown that BCS superconductors have peaks in the 
second derivative of the absorption spectrum {A = 1 — i?, 
where R is the reflectance)EirE3 at the frequencies of lon- 
gitudinal and transverse phonons, thus confirming the 
mechanism of electron-phonon interaction in these ma- 
terials. For this test one has to be able to measure the 
second derivative of the reflectance in the region of the 
relevant excitations. To get the good signal-to-noise ratio 
necessary for such experiments large single crystals are 
needed. 

Several studies of the opticat-jpitdjictivity of the rbpro- 
carbides have been published S-ESEa Widder et alt3 re- 
ported on reflectivity and electron loss spectroscopy on 
ceramic samples of LuNi2B2C at room temperature for 
frequencies up to 50 eV. They find an overall metallic 
response with a plasma frequency of 4.25 eV and an 
electron-phonon coupling parameter Xtr = 1.2 . Bom- 
meli et al. measured both LuNi2B2C and YNiaJB2C over 



a wide range of temperatures and frequenciesO report- 
ing a strong absorption in the pure LuNi2B2C material 
with an onset of 100 cm^^, well below any significant 
phonon density of states. They also claim to see a dirty 
limit superconducting gap in the undoned material while, 
according to transport measurements Jj the undoped ma- 
terials should be in the clean limit. In a recent paper 
Kim et a/.|-show measurements at room temperature in 
LuNi2B2CE3 supplementing infrared reflectance with el- 
lipsometry measurements at higher frequency. They also 
studied the effect of annealing on their polished samples. 
They reported a screened plasma frequency of 3.76 eV for 
an annealed sample and 3.01 eV for an unannealed sam- 
ple. Both samples were polished. They find a Drude-like 
frequency-independent scattering rate below 0.3 eV but 
a gradual rise in scattering above this frequency signaling 
non-Drude behaviour in the midinfrared. 

Cheon et al. reported on dc transpO|Ct in a series of 
cobalt doped LuNi2B2C single crystals.El They found a 
dc resistivity that was linear in temperature with an ia^ 
tercept on the temperature axis at 35 K at T = 0£)/lO.E3 
The slope of the resistivity increased only slightly with 
cobalt doping, providing evidence that any changes in the 
density of states at the Fermi surface with cobalt doping 
were small. We decided to investigate the issue of the 
superconducting gap in clean vs. dirty samples in the 
borocarbides using samples from the Ames group, the 
same source as those used in the Bomelli et al. work. 



II. EXPERIMENT 

We used flux grown single crystals in this work char- 
acterized extjeftsjvely by dc resistivity and magnetic 
susceptibility.BEa The reflectance measurements were 
done on natural growth faces normal to the c-axis. To 
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FIG. 1: The reflectance of LuNi2B2C at two temperatures. 
The data above 5 eV is from Widder et a/E3 



correct for irregularities of the surface, which were mainly 
vertical steps, the sample reflectance was referenced to 
spectra where the sample was coated with a gold layer, 
evaporated in situBa The as-grown crystals have tiny 
droplets of Ni2B on them. The visible ones occupy a 
small fraction of the surface area and our gold evapora- 
tion technique tends to suppress curved surfaces in favour 
of flat ones. Metallic droplets that are much smaller than 
the wavelength can give rise to an absorption band in the 
region of a;p/\/3. 

Most of the measurements were carried out in a cold- 
finger flow cryostat with a minimum sample temperature 
of 10 K. A home made rapid-scan interferometer was used 
in the far-infrared with He-cooled bolometer detectors 
below 700 cmr^ and an MCT detector up to 8000 cm~^ 
High frequency measurements were done with a grating 
spectrometer, and in the superconducting state, a polar- 
izing interferometer and an immersion dewar with a He-3 
bolometer were used between 5 and 60 cm~^. 

Figure |^ shows the reflectance of a crystal of pure 
LuNi2B2C over a wide range of frequencies at ambient 
temperature and at 14 K. Above 5 eV we have merged 
our data with the EELS data of Widder et aZ.E3 The data 
showed little temperature dependence at 8000 cm~^ and 
we made only room temperature measurements above 
this frequency. At low frequency the reflectance in the 
different spectral regions, measured with different spec- 
trometers and detectors, generally agreed to within 0.5 
% in the region of overlap, whereas at higher frequencies 
the mismatch could be as much as several percent. 

Figure |^ shows the low frequency region on an ex- 
panded scale at three temperatures. The solid lines refer 
to the pure sample of LuNi2B2C and the dashed curves 
to the sample with 9 % Co doping. The fine structure on 
the curves is noise. 

A comparison of our measured refiectance with data 
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FIG. 2: Low-frequency reflectance of LuNi2B2C at three tem- 
peratures. The sohd curves are for the undoped sample and 
the dashed ones for a sample with 9 % cobalt doping. The 
fine structure is noise. 



from other groups shows in general reasonably good 
agreement with most of the previous work but with some 
notable discrepancies. First, our reflectances are in gen- 
eral higher than what has been reported by previous in- 
vestigators. Part of this difference we attribute to our use 
of as-grown single crystal surfaces since polishing will in- 
troduce surface damage leading to additional absorption 
particularly at higher frequencies. Kim et al. find that by 
annealing the polished crystals a substantial increase in 
reflectance can be obtained. cJ We have found, in the past, 
that to obtain good reflectance spectra of metals in the. 
ultraviolet it is necessary to electropolish the samples.ca 
The reflectances of our as-grown single crystal surfaces 
are much higher than the polished samples leading to 
conductivities in the interband region that are as much 
as three times higher than those of the polished samples. 

At 800 cm~^ at room temperature, we measure a re- 
flectance of 94 % whereas both Kim et al. and Bommeli 
et al. report ^91 % in polished samples. We fail to see 
the strong absorption seen by Bommeli et al. at low tem- 
perature with an onset at 50 cm^^ where the reflectance 
falls rapidly from 99 % at 50 to 93 % at 150 cm"!. In all 
our samples, doped and undoped, the low temperature 
reflectance below 700 cm~^ was above 98 %, as expected 
for a good metal. 

The overall high reflectance of these single crystals 
places demands an an accurate calibration of the abso- 
lute value of the reflectance. We calibrated our system by 
using polished stainless steel and a lead indium alloy as 
primary standards. We started by measuring the dc resis- 
tivity of these reference samples and used a Drude model 
to calculate the absolute value of the reflectance. These 
calibrations were used to determine the reflectance of our 
evaporated gold coating which was our secondary stan- 
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FIG. 3: The temperature dependence of the reflectance at 
two frequencies. The circles are data for the pure sample, 
the squares for the cobalt-doped sample. The dashed lines 
are least squares fits to the experimental points. The small 
points are for a theoretical model based on the interaction 
with a mode at 170 cm~^ with parameters fit to dc transport 
and an assumed plasma frequency of 3.3 eV (top points pure 
and bottom doped samples). 



dard. Runs were rejected where the reflectance changed 
by more than 0.5 % between the beginning and the end of 
the run. As another check we compared our reflectances 
in the low frequency region with the Drude model based 
on dc resistivity measurements of crystals prepared in the 
same way as shown in Fig. ^ and Fig. ^. 

Figure ^ gives a picture of the reproducibility of our 
measurements. We have plotted the reflectance at 350 
cm^^ and at 700 cm~^ as a function of temperature for 
our two samples. In this spectral region lut > 1 and 
R — 1 — 2/ (ujpT). As the discussion below shows, we find 
a plasma frequency of LUp = 3.3 eV. With this value we 
can relate our reflectance measurements directly to dc re- 
sistivity {p = An/ (WpT)). The predicted reflectances from 
the dc resistivity are shown as small symbols in Fig. ^jfor 
the two samples. We see generally excellent agreement in 
the temperature dependence as well as the difference in 
the two sample curves due to the doping induced elastic 
scattering. The dashed curves are a least squares fit of a 
straight line to the experiments. The root mean square 
deviation from the straight line is less than 0.15 % for 
all the curves. However the high frequency curves fall 
systematically below the dc transport prediction. This is 
due to the rising 0.15 eV mid-infrared band that has a 
threshold in this spectral region. 

To obtain the optical conductivity we performed 
Kramers-Kronig analysis of the refiectance spectra. For 
this it is necessary to extend the data at low and high fre- 
quencies. At low frequencies we assumed a Drude model 
where the parameters were fit to the lowest measured fre- 
quency, 60 cm~^ for the undoped sample and 20 cm~^ 
for the cobalt doped sample. At high frequency, £mm 5 
eV to 35 eV, we used the data of Widder et alr3 and 
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FIG. 4: The optical conductivity of pure LuNi2B2C at three 
temperatures. The dashed lines show estimated error limits 
based on a 0.5 % uncertainty of the 100 % line. The circular 
symbol on the vertical axis is the dc conductivity at room 
temperature. The decreasing conductivity as the temperature 
is lowered is due to the narrowing of the Drude peak. 



FIG. 5: The conductivity of the 9 % cobalt-doped sample 
(dashed lines) and the undoped sample, (solid lines). The 
solid ball on the vertical axis is the dc conductivity. While 
there is little difference at room temperature, the Drude peak 
of the doped sample is substantially broader at low tempera- 
ture. 



beyond that a power law R (x uj ^■'^ to uj — 100 eV and 
R cx uj-'^ above 100 eV. 

Figure ^ shows the optical conductivity for the un- 
doped sample at 300, 200 and 14 K for the whole mea- 
sured range. The dashed curves are estimated errors 
based on a 0.5 % error in the absolute value of reflectance. 
The solid point is the dc resistivity measurement at room 
temperature of Cheon et alu Our optical data match the 
dc data within our experimental error. The error of the 
dc data is ±10% and mainly due to uncertainties in con- 
tact geometry. As the temperature is lowered the con- 
ductivity drops in the midinfrared and rises in the far 
infrared. We are unable to show far infrared conductivi- 
ties below 200 cm~^ because the reflectance becomes too 
high for accurate Kramers-Kronig analysis. 

Figure || shows the conductivity of the doped sample 
at two temperatures, along with the pure sample, shown 
with solid lines. It can be seen clearly that while the 
doping has little effect at room temperature, at low tem- 
perature the Drude peak is broader as expected from the 
additional scattering. 

The measurements with the polarizing spectrometer in 
the superconducting state for the 9 % doped sample are 
shown in Fig. ^. To increase the sensitivity we only show 
a reflectance ratio between the superconducting state and 
normal state shown as the solid curve in the figure. The 
dirty limit BCS superconducting gap signature can be 
seen clearly as a drop in the ratio at 22 cm^^. The dashed 
curve shows the calculated reflectance ratio based on dc 
conductivity of the 9 % doped sample. We estimate a 
gap value 2A = 24 ± 2 cm~^. 

The disagreement in our data of the peak amplitude 
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FIG. 6: The reflectance in the superconducting state at 1.3 
K divided by the normal state reflectance at 13.5 K. The 
superconducting gap can be seen as an onset of absorption in 
the normal state. A BCS theoretical curve for this sample is 
shown as the curve with the long dashes. The dotted curve is 
a calculation for the undoped sample. The undoped sample 
is in the clean limit and the gap signature is predicted to be 
weak. 



in the observed reflectance ratio with a BCS calculation 
is 0.5 %. This difference is larger than we expect from 
such a measurement where the sample is not moved and 
the two measurements are done within minutes of one 
another. We have measured a lead indium alloy in the 
same geometry without observing a discrepancy between 
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the dc transport prediction and the reflectance ratio. The 
calculation assumes the superconductor to be a perfect 
reflector in the gap region. Any residual absorption in the 
superconducting state due to strong gap anisotropy or a 
second phase on the sample surface would reduce the ab- 
sorption below the calculated value and have the opposite 
effect to what we observe. An overall decrease in detector 
sensitivity with increasing temperature may be partially 
responsible for the discrepancy but it would have the ef- 
fect of uniformly shifting the ratio upward causing the 
data above the gap value to exceed unity. We do not 
observe such an effect here although we cannot rule out 
a 0.2 % contribution from this cause. Since we have used 
dc resistivity to estimate the strength of the normal state 
absorption any discrepancy between the resistivity of our 
sample and the bulk measurement could explain the dis- 
crepancy. However the agreement at higher frequencies 
shown in Fig. ^ rules out any large discrepancies. 

The model predicts a sharp change of slope at 2A 
whereas our experiments show what might be two steps, 
one at 22 cm-1 and another at 26 cm-1, which way be due 
to gap anisotropy as proposed by Yang et alB based on 
Raman spectroscopy done in different scattering geome- 
tries. However scattering due to the Co dopants would 
be expected to remove any such anisotropy. 

The curve with the long dashes shows the calculated 
spectrum for the pure sample. The expected gap signa- 
ture is weak since the sample approaches the clean limit 
with l/r = 17 cm~^, determined from the dc conductiv- 
ity with ujp — 3.3 eV and 2A = 46 cm~^ from Raman 
spectroscopy. We measured the pure sample reflectance 
ratio but failed to see any sign of a superconducting gap 
in the reflectance ratio of our samples within the esti- 
mated noise level of ± 0.002 in the reflectance ratio. 



III. RESULTS 

There are two popular ways of analyzing the conduc- 
tivity if it does not follow the standard Drude form. The 
first, the two component method, is to fit the conduc- 
tivity to a Drude peak centered at zero frequency and 
add a second component consisting of a set of oscilla- 
tors in the mid-infrared to allow for any additional non- 
Drude absorption. The second method, the exteftded 
Drude model, introduced by Allen and Mikkelsen,EZl as- 
sumes that the Drude form applies throughout the in- 
frared but that the Drude scattering rate increases with 
frequency, making, in effect, the Drude peak broader as 
the frequency increases. The first method accounts well 
for any parallel channels of conductivity such as optic 
phonons or low-lying interband transitions. The second 
method is needed when the interaction of the free carriers 
is strong and a Holstein sideband is formed at low tem- 
peratures. At high temperatures this sideband merges 
with the Drude absorption forming an approximate over- 
all Drude-like band. In the absence of a detailed theory 
of the optical conductivity both methods have value in 
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FIG. 7: Oscillator fits to the reflectance (solid curve). In 
addition to the Drude term at low frequency, three oscillators 
are needed to describe the conductivity at higher frequencies. 



the analysis of the data if one keeps open the possibil- 
ity that some of the mid-infrared oscillators have no real 
physical meaning or that the increase in the scattering 
rate seen in the one-component analysis may be due to 
interband processes and not to self energy effects on the 
free carriers. 

The two component picture is illustrated in Fig. |^ 
where we show an oscillator fit to the conductivity for the 
undoped sample at room temperature. For the undoped 
sample at 300 K the following oscillators were used, (cen- 
ter frequency, width, and plasma frequency, in cm~^): 0, 
500, 25200; 1200, 1200, 15000; 8410, 20000, 60000; 90000, 
150000, 180000. It is clear that the simple Drude form 
does not fit at high frequency and additional Lorentz 
oscillators are needed to fit the data. The two high fre- 
quency oscillators can be identified as due to interband 
transitionfO since band structure calculations indicate a 
gap in the density of states just below the Fermi level 
rising to a strong peak about 2 eV below the Fermi level. 
Combining this hole density of states with the sharp peak 
at the Fermi level would predict an optical conductiv- 
ity that rises gradually from low frequencies to form the 
plateau we see in the 1 to 10 eV region. It is difficult, 
in the absence of a detailed joint density of states, to 
say if there is a low frequency threshold to the optical 
conductivity but since the dispersion curves show several 
band crossings very close to the Fermi level, one cannot 
rule out an interband conductivity that extends to zero 
frequency. 

The midinfrared oscillator at 0.15 eV is more difficult 
to assign. Widder et alE3 have assigned the spectral 
weight in this region to the frequency dependent scat- 
tering rate of phonons, but as we will show below, it is 
difficult to get a consistent picture of the data in this 
spectral region with a phonon spectrum that is confined 
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to the known phonon density of states. Fig. ^ (to be dis- 
cussed below) shows that there is an absorption threshold 
at 700 cm^^ at low temperature which suggesting that 
interband transitions may set in at this frequency. How- 
ever the onset frequency of the threshold is strongly tem- 
perature dependent, moving to 100 cm^^ at 300 K. This 
temperature dependence is difficult to explain in terms 
of a simple Fermi function expected for interband tran- 
sitions where a broadening of only fc^T w 200 cm~^ is 
expected. 

The spectral weight under the 0.15 eV oscillator could 
also be the result of a multiphonon side-band of the 
Drude peak. From its spectral weight we estimate an ad- 
ditional mass-enhancement factor of 1.35 and a total free 
carrier plasma frequency of 3.64 as opposed to 3.12 eV if 
the 0.15 eV oscillator is assigned to interband transitions. 
A complete Eliashberg calculation would be needed to see 
if this model will give a consistent picture of the gap ratio, 
the coupling constant A and the transition temperature 
Tc. 

There are several other ways of estimating the Drude 
weight from our data. The Drude component in Fig. |^ 
has a plasma frequency of 3.14 eV for the undoped sam- 
ple. If we scale it to agree with the dc resistivity we 
get 3.6 eV. Another method uses the temperature de- 
pendence of the reflectance shown in Fig. |^ and com- 
pares it to the temperature dependence of the dc resis- 
tivity. The derivative of the dc resistivity is given by 
dp/dT = ^d-y/dT where 7 is the scattering rate. On 

the other hand the reflectance in the relaxation region, 
where 7 < w is given by 1 — i? = 27/wp, where R is 
the reflectance. This equation can be differentiated with 
respect to T and the relaxation rate derivative d-y/dT 
can be eliminated between these two equations to give 
ujp — 27r(l — R)' / p' where the primes denote tempera- 
ture derivatives. This method is not subject to errors 
in absolute reflectance as it uses the temperature varia- 
tion as input. Using this method with the fitted curves 
in Jig. 1^. along with the resistivity data of Cheon et 
alM we obtain an average plasma frequency of 3.4 ± 0.15 
eV directly from the reflectance without using data from 
Kramers-Kronig analysis. As an average of the two meth- 
ods we adopt the plasma frequency of our pure sample 
of 3.3 ± 0.2 eV. 

Our results are lower than the 4.25 eV value of Widder 
et al. who included the 0.15 eV peak in the Drude spec- 
tral weight. Bommeli et al. find a two-component low- 
frequency spectrum with a Drude weight of 0.73 eV and 
an additional mid infrared band with a total oscillator 
strength of 4 eV by including the broad absorption at 500 
cm~^. We have found no evidence for this low frequency 
component. Kim et al. using polished and annealed sam- 
ples find LUp = 3.76 eV for the screened plasma frequency 
determined from the zero crossing of ei. The zero cross- 
ing of ei, obtained from Kramers-Kronig analysis of our 
reflectance occurs at the same frequency. However to get 
the unscreened plasma frequency Up a correction has to 
be made for the very high frequency dielectric constant 
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FIG. 8: Frequency dependent scattering rate, top panel, and 
the effective mass, lower panel, obtained from the extended 
Drude formula, assuming a plasma frequency of 3.3 eV. The 
dashed line is for the cobalt-doped sample. The symbols on 
the axis are dc values of the scattering rate, corresponding to 
the optical conductivity curves. 



which is not known accurately. 

The frequency-dependent scattering rate formalism 
starts with the extended Drude formula: 



aiiu,T) = 



1 



4tT 1/t{uJ, T) - + Xtr{iO, T)] ' 



where 1/t{uj,T) is the frequency dependent scattering 
rate and 1 -I- Xtr{i-^,T) — m* /m is the mass renormaliza- 
tion factor. These can be calculated from the real and 
imaginary parts of the optical conductivity: 



and 



1 + XtrH 



1 



The frequency dependent scattering rate is shown in 
Fig. H along with the mass renormalization. Both are 
based on an assumed plasma frequency of 3.3 eV for the 
free carriers and a high frequency dielectric constant of 
1.0. The choice of plasma frequency affects the overall 
scale of the curves. The choice of high frequency di- 
electric constant has little effect in this frequency region 
where ei and £2 are large. We see a monotonic increase in 
scattering rate with frequency. The solid curves are for 
the pure sample and the dashed one for the doped one. 
The symbols on the vertical axis are the values from the 
dc conductivity assuming a plasma frequency of 3.3 eV. 

The lower panel of Fig. ^ shows m*/m= 1 + Atrforthe 
pure sample at low temperature, which is, within a factor 
of cj, the imaginary part of l/r. One can read off the cal- 
culated electron-phonon coupling constant for dc trans- 
port, the zero-frequency value, Atr.(0) to be 0.15. The 
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value obtained more directly from dc transport is 0.34, 
based on a plasma frequency of 3.3 eV and the assump- 
tion that the transport curves have reached their high 
frequency linear limit at 300 K when \/t — 2TTkBT. A 
strong contribution from high frequency phonons would 
raise the dc value of dp/dT and lead to an even larger 
At,-. We do not feel that discrepancy between the trans- 
port value and the optical value is significant since the 
curves in Fig. ^ have large uncertainties due to the high 
reflectance of the undoped sample at low frequencies. 

It is important to stress the difference between the 
curves iii— Fig. ^ and the corresponding data for the 
cupratesEj First, the overall scale of scattering in the 
cuprates is an order of magnitude higher to the point 
where the scattering rate is larger than the frequency, 
leading to a breakdown of Fermi liquid theory which 
assumes that l/r <^ lo. Second, while the variation 
of scattering is linear with frequency in both materials, 
in the cuprates the curves extrapolate back to a finite 
positive intercept on the l/r axis that scales as T and 
goes to zero at "F-i^ 0. This is characteristic of quan- 
tum fluctuationsEil For the borocarbides the intercept at 



T = is on the negative l/r side resulting in a finite 
intercept on the frequency axis, which is of the order of 
the Debye frequency, as expected for a Bloch-Griineisen 
conductivity. The Debye frequency of LuNi2B2C is 250 
cm^^ from a Debye temperature of 360 K measured by 
neutron spectroscopjilj in good agreement with the point 
of intersection of the low temperature reflectance curve 
with the R = \ axis. 

A fit of a theoretical model directly to the reflec- 
tivity avoids the propagation of errors associated with 
Kramers-Kronig analysis and the extraction of the pa- 
rameters of the extended Drude model. This is par- 
ticularly relevant in a good metallic system where the 
reflectivity is close to unity. The procedure is to start 
with a niode] al^{Vl)F{Q) of the electron phonon spectral 
functionHL^ weighted by the amplitude for large angle 
scattering and then calculate the optical properties and 
compare them to experiments. 

To calculate the reflectance spectra we use the ex- 
tended Drude fjetmalism with the scattering rate given 
by Shulga et alEj: 



-{lo, T) = -j^ dnal{n)F{n)[2u:coi\i{ — ) ~ [u + f])coth(^^) + {uj - r!)coth(^^) 



where T is the temperature, measured in frequency units. 
Following Gonnelli et alx^ we place a restriction on the 
spectral functions demanding that they fit the slope and 
absolute magnitude of the dc resistivity of LuNi2B2C. 
We use a plasma frequency of 3.3 eV in our models. 

The simplest model is a single Einstein oscillator. We 
find that an oscillator at 170 cm^^ with \tr — 0.33 gives a 
fairly good value of the room temperature resistivity and 
the temperature derivative of the resistivity (35 ^ficm 
and 0.135 /ificm/K, vs. the experimental values of 36 
/ificm and 0.127 /xr2cm/K). The bold curves in Fig. ^ 
are the calculated spectra based on this model. 

Another model uses the phonon density of states as 
determined by neutron scattering. The bare spectrum 
gives a rather poor fit to the dc transport and Gonnelli 
et al. found that they had to enhance the low frequency 
modes (below 38 meV) by 0.7 and the high frequency 
ones by 0.25, a ratio of 2.8 in the enhancement of the 
low frequency spectrum. They used a \tr of 0.53 based 
on a plasma frequency of 4.25 from the Widder et alxn 
optical measurements. 

Because we find a lower plasma frequency of 3.3 eV 
we find that somewhat different parameters than those 
used by Gonnelli et al. are necessary to fit the resis- 
tivity. Taking the Gonnelli spectrum as a reference, we 
have multiplied the low frequency part by 0.23, the upper 
band above 38 meV by 0.48, a ratio of 2.1 in the enhance- 



ment of the low frequency spectrum. These parameters 
were adjusted to give reasonable agreement with the tem- 
perature dependence of the dc resistivity. The resulting 
spectrum is similar to what is shown for the Einstein 
oscillator in Fig. ^. 

The fit of both models to our reflectance is good at 
low frequency and well within experimental uncertain- 
ties, but there are serious deviations at high frequency 
which become more marked with increasing temperature. 
These deviations exceed our experimental uncertainty. 

It is not possible to remove the high frequency discrep- 
ancy by manipulating the coupling to the various phonon 
frequencies. To obtain the strong scattering seen at high 
frequencies it is necessary to add oscillators above the 
phonon spectrum with a continuous distribution of fre- 
quencies, a spectrum of the Marginal |Eermi liquid type 
as seen in the cuprate superconductors£j However such a 
model cannot be justified in this case since the scattering 
rate curves shown in Fig. ^ intersect above unity, signi- 
fying the coupling to finite frequency bosons and not to 
quantum critical fiuctuations where the-intercept occurs 
at unit reflectance at low temperature.Eil 

We think a more likely scenario is one where the thresh- 
old of absorption seen in Fig. || is a manifestation of the 
onset of either multi-phonon processes or interband tran- 
sitions as shown by the 0.15 eV band in the conductivity 
spectra, the same phenomenon that was responsible for 
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FIG. 9: Calculated reflectance based on the electron-phonon 
interaction with a coupling constant Xtr obtained from the dc 
resistivity assuming an Einstein spectrum of a single phonon 
at 170 cm~^. The agreement with the measured reflectance 
(thin curves) is good at low frequency but strong additional 
absorption sets in at 700 cm~^ at 14 K 

the mid-infrared oscillator shown in Fig. ^. There is 
also the possibility of a contribution to the mid-infrared 
absorption by the Ni2B droplets. 

IV. DISCUSSION 

We first address the question of the plasma frequency 
and its variation with doping. Based on several methods 
of analysis we find a Drude weight at room tempera- 
ture of 3.3 eV for the undoped sample. Other investiga- 
tors have reported somewhat higher values as a result of 
the inclusion of some midinfrared spectral weighJ_in the 
Drude peak gpfing a value closer to 4.25 eV.EjO Pick- 
ett and Singhli3 calculate a plasma frequency -fipm band 
structure to be 5.1 eV. However Michor et al!3 suggest, 
from an analysis of specific heat data, that this is an over 
estimate. 

The question of the doping dependence of the Drude 
weight is more difficult. It is clear from most of our mea- 
surements that the changes in the spectra with Co doping 
are relatively slight and the methods used in the previ- 
ous paragraphs are subject to errors of the order of 10% 
and not accurate enough for this task. To help reduce 
systematic errors we use the integrated spectral weight 
as a measure of Drude weight. Fig. ^ shows the quan- 
tity Nef f which is the number of electrons contributing 
to the conductivity in the unit cell of volume Vceii ■ 

^eff{^) = — / cr(w )(Lj 

We show Neff at room temperature, for the pure sample, 
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FIG. 10: Partial conductivity sum rule for the pure and the 
doped samples. The small enhancement of conductivity in the 
pure sample would imply a small enhancement of the density 
of states at the Fermi surface. The small difference is not 
significant within our margin of error. 



solid curve, and the 9 % doped sample, dashed curve. 
The pure curve is higher by about 7 % in the region of 
the Drude peak. This difference is not significant in view 
of our error of 0.5 % in reflectance which propagates to 
give an uncertainty of 7 % for the spectral weight. The 
density of states at the Fermi surface is proportional to 
the square of the plasma frequency and appears to be 
reduced by 7 % by the cobalt doping if we believe the 
difference data of Fig. |l^. The overall spectral weight 
reaches one electron per unit cell at the frequency where 
the interband contribution takes over from the Drude 
conductivity. 

We can compare our change in Drude spectral weight 
with cobalt doping with what might be expected from the 
dc resistivity. From the data of Cheon et alu the slope 
of the resistivity changes from 0.127 to 0.140 ^fJcm/K in 
going from the pure sample to the 9 % cobalt doped one. 
If this effect is from the change in the Drude plasma 
frequency it would correspond to a change in uj^ of 10 
%, somewhat higher than the 7 % that we get from the 
change in the infrared Drude weight. However the dc 
measurements are subject to an error of ±10 % due to 
uncertainties in contact geometry and the difference is 
not significant. It should be noted that the dc resistivity 
slope is also proportional to the coupling constant Xtr 
which could change with doping. It appears that neither 
method, the infrared or the dc transport is able to put an 
accurate limit on the change in Drude weight with cobalt 
doping. All that can be said is that it appears to be less 
than 10 % from both measurements. 

Our value of 2A/fcBTc = 3.9 ± 0.4 for the doped 
sample places it in the regime of weak to moderately 
strong coupling superconductivity. Raman scattering on 
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a pure sample of this material gives a strong coupling 
gap ratio of 4.1,0 as do thermodynanaic data, such as the 
magnitude of the specific heat jumpJ13 However, tunnel- 
ing spectroscopy point to a weak coupling value of the 
ratioEil A weak coupling value of the gap parameter is 
also consistent with our electron phonon coupling param- 
eter Xtr = 0.33.eI 

In summary, we have presented reflectance data on 
as-grown surfaces of LuNi2B2C which differ in several 
aspects from those measured by previous investigators 
along with an analysis based on free electron theory. First 
our measurements confirm that on the whole this mate- 
rial has the electrodynamic properties of a good metal. 
We are able to fit the reflectance data at low frequency 
to models of electron phonon interaction with parame- 
ters derived from dc transport on crystals from the same 
source. At higher frequency additional absorption is 
observed which we attribute to multi-phonon processes 
or interband transitions. We find a Drude plasma fre- 
quency of 3.3 eV and an electron-phonon coupling con- 
stant Xtr = 0.33. This low value of Xtr is consistent 
with our superconducting gap ratio of 2A/kBTc = 3.9 
which we measure from reflectance in the superconduct- 
ing state in the cobalt doped sample which is in the dirty 
limit. We do not see any change to the plasma frequency 
with cobalt doping to a level of ± 10 % but do observe an 
amount of increased scattering which is consistent with 
Mattheissen's rule. 

There remain several open questions. First, what is 



the nature of the mid-infrared absorption band at 0.15 
eV? Unlike the cuprates where a large portion of the 
absorption is associated with a Marginal Fermi Liquid 
spectrum of electronic coupling to the charge carriers 
with an intercept at zero frequency, the intercept here 
is Bloch-Griineisen-like at finite frequency. Second, the 
accuracy of our experiments was insufficient to yield the 
detail in scattering rate spectrum necessary to reveal the 
phonon spectrum. With more accurate infrared data the 
a^F{fl) spectrum of the electron phonaa interaction can 
be mapped out as was done for lead.E3 Also, it would 
be important to make accurate resistivity measurements 
at high temperature to look for evidence of coupling to 
high frequency modes. To see evidence for phonon ef- 
fects by infrared spectroscopy, the current single bounce 
reflection measurements do not have enough resolution. 
Cavity measurements on large enough single crystals at 
low temperature and high frequency may have enough 
sensitivity for this. Also, a detailed strong coupling cal- 
culation may be able to throw light on the possibility 
that the 0.15 eV band may be due to strong coupling 
multi phonon processes. 
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